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A STUDY OF THE MOTION AND AERODYNAMIC HEATING 
OF MESSIIES ENTERING THE EARIH'S ATMOSPHERE 
AT HIGH SUPERSONIC SPEEDS 
By H. Julian Allen and A. J. Eggers, Jr. 


SUMMARY 


A simplified analysis is made of the velocity and deceleration 
history of missiles entering the earth's atmosphere at high supersonic 
speeds. It is found that, in general, the gravity force is negligible 
con 5 )ared to the aerodynamic drag force and, hence, that the trajectozy 
is essentially a straight line. A constant drag coefficient and an 
exponential variation of density with altitude are assumed and general- 
ized curves for the variation of missile speed and deceleration with 
„ altitude are obtained. A curiotis finding is that the maxi mum decelera- 

tion is independent of physical characteristics of a missile (e.g., mass, 
size, and drag coefficient) and is determined only by entry speed and 
*• flight-path angle, provided this deceleration occurs before impact. 

This provision is satisfied by missiles presently of more usual interest. 

The results of the motion analysis are employed to determine means 
available to the designer for minimizing aerodynamic heating. Emphasis 
is placed upon the convective -heating problem including not only the 
total heat transfer but also the maximum average and local rates of 
heat transfer per unit eirea. It is found that if a missile is so heavy 
as to be retarded only slightly by aerodynamic drag, irrespective of 
the magnitude of the drag force, then convective heating is minimized 
by minimizing the total shear force acting on the body. This condition 
is achieved by employing shapes with a low pressure drag. On the other 
hand, if a missile is so light as to be decelerated to relatively low 
speeds, even if acted upon by low drag forces, then convective heating 
is minimized by employing shapes with a high pressure drag, thereby 
maximizing the amount of heat delivered to the atmosphere and minimizing 
the amount delivered to the body in the deceleration process. Blunt 
» shapes appear superior to slender shapes from the standpoint of having 

lower maximum convective heat-transfer rates in the region of the nose. 
The maximum average heat-transfer rate per unit area can be reduced by 
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employing either slender or hlunt shapes rather than shapes of inter- 
mediate slenderness. CJenerally, the blnnt shape vith high pressure 
drag wonld appear to offer considerable promise of minimizing the heat 
transfer to missiles of the sizes, weights, and speeds presently of 
interest. 


INTRODUCTION 


In the design of long-range rocket missiles of the ballistic type, 
one of the most difficvilt phases of fli^t the designer must cope with 
is the re-entry into the earth's atmosphere, wherein the aerodynamic 
heating associated with the high flight speeds of such missiles is 
intense. The air temperat\ire in the boundary layer may reach values in 
the tens of thousands of degrees Fahrenheit which, combined with the 
high surface shear, promotes very great convective heat transfer to the 
surface. Heat-absorbent material must therefore be provided to prevent 
destruction of the essential elements of the missile. It is a charac- 
teristic of long-range rockets that for every pound of material which 
is carried to ”bum-out,” many povinds of fuel are reqvdLred in the 
booster to obtain the flight range. It is clear, therefore, that the 
amount of material added to protect the warhead from excessive aero- 
dynamic heating must be minimized in order to keep the take-off weight 
to a practicable value. The importance of reducing the heat transferred 
to the missile to the least amount is thus evident. 

For missiles designed to absorb the heat within the solid surface 
of the missile shell, a factor which may be importeuit, in addition to 
the total amount of heat transferred, is the rate at which it is trans- 
ferred since there is a maximum rate at whiqh the surface material can 
satfely conduct the heat within itself. An excessively high time rate 
of heat input may promote such large tenq^erature differences as to 
caiise spalling of the surface, and thus result in loss of valuable heat- 
absorbent material, cac even structural failure as a result of stresses 
induced by the tenr^rature gradients. 

For missiles designed to absorb the heat with liquid coolants 
(e.g., by "sweat cooling” where the surface heat-transfer rate is high, 
or by circvilating liquid coolants within the shell where the surface 
heat-transfer rate is lower) , the time rate of heat transfer is simi- 
larly of interest since it deterxaines the required liquid ptnnping rate. 

These heating problems, of course, have been given considerable 
study in connection with the design of particular missiles, but these 
studies are very detailed in scope. There has been need for a general- 
ized heating analysis intended to show in the broad sense the means 
available for minimizing the heating problems. Wagner, reference 1, 
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made a step toward satisfying this need "by developing a laudably simple 
motion analysis. This analysis was not generalized, however, since it 
was his purpose to study the motion and heating of a particular missile. 

It is the piorpose of this report to simplify and generalize the 
analysis of the heating problem in order that the salient feat\ires of 
this problem will be made clear so that successful solutions of the 
problem will suggest themselves . 

A motion analysis, having the basic character of Wagner’s approach, 
precedes the heating analysis. The generalized results of this analysis 
are of considerable interest in themselves and, accordingly, are treated 
in detail. 


ANALYSIS 

Motion of the Body 


Consider a body of mass m entering the atmosphere from great 
height. If, at any altitude y, the speed is V and the angle of 
approach is 0 to the horizontal ( see sketch) , the parametric equa- 
tions of motion can be written^ 


d5y 

dt^ 


dt^ 


Ct^PV^A 

- g + sin 0 


CjjPV^A 

2m 


2m 


cos 0 



^Properly, the analysis should consider those effects resulting from the 
fact that the earth is a rotating sphere, but since the altitude range 
for which drag effects are Important is less than 1 percent of the 
radius of the earth, the rectilinear treatment given in this analysis 
is permissible. 
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vhere 

Cj) drsig coefficient, dimensionless 
V speed, ft /sec 

A reference area for drag evaluation, sq, ft 
m mass of the tody, slugs 

P mass density of the air, slvigs/ft® 

g acceleration of gravity, ft/sec^ 

x,y horizontal and vertical distance from the point of iii 5 >act with 
the earth, ft 

9 angle between the flight path and the horizontal, deg 

(See Appendix A for con^>lete list of symbols.) 

In general, the drag coefficient varies with Mach ntanber and 
Reynolds number, while the density and, to a very minor extent, the 
acceleration of gravity vary with eiltitude. Hence it is clear that 
exact solution of these equations is formidable. Let us first, then, 
consider the following simplified case: 

1. The body descends vertically. 

2. The drag coefficient is constant.^ 

3. The acceleration of gravity is constant.® 

i)-. The density as a function of eiltitude is given by the relation 

0 = ( 2 ) 

where Pq and p are constants^ 

^As is well known, this assumption is generally of good accuracy at the 
high Mach numbers under consideration, at least as long as the total 
drag is largely pressure drag. 

®The acceleration of gravity decreases by only 1 percent for every 
100,000-foot increase in altitude. 

‘^‘This relation is consistent with the assumption of an isothermal 
atmosphere. 
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Equations (l) then reduce to the single equation 


dt^ dt 2m 


( 3 ) 


Noting that 


we let 


^ _ V — 

dt dy 


Z = V" 


and equation (3) hecomes the linear differential equation 


3£_^e-«i^Z + 2s = 0 

dy HI 


(1^) 


which has the weH-knovm solution 




‘^D^o-^ -Py, 
e ^"^dyi 


m 


Z = e 


- 2g / e 




-/ 


CpfoA -py 


m 


e '■^^dy 


dy + const. 


Performing the Integrations, we obtain as the solution of this relation 


-for 


Z = V"" = e 


pm 


\n 


w 

£s V 
p L 


n=i 


n]n 


- 2gy + const. 


(5) 


so that the deceleration becomes, in terms of gravity acceleration. 


dV 




g 2mg 


9pPcA - 

pm. 


iSy 


w 


n=i 


n[n 


- 2gy+ const. 


-1 

(6) 
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As an example, consider the vertical descent of a solid iron sphere 
having a diameter of 1 foot. For a sphere the drag coefficient may he 
taken as unity, based on the frontal area for an Mach manhers greater 
than about l.^l-. In eq.Taation (2), irtiich describes the variation of 
density with altitude, the constats should clearly be so chosen as to 
give accurate values t»f the density over the range of altitudes for 
which the deceleration is large. It is seen in figure 1 that for 


= 0.003l{- slugs /ft° 

and 

P = — ft"^ 

22,000 

which yields 

y 

P = 0.003^<- e ^^^Tooo 

the calculated density is in good agreement with the NACA standard 
atmosphere val\ies obtained from references 2 and 3 for the altitude 
range from 20,000 to l80,000 feet. These relations have been tised in 
calculating the velocity and deceleration of the sphere for various 
altitudes, assuming vertical entrance velocities of 10,000, 20,000, 
and 30,000 feet per second at ItO miles altitude which, for these cases, 
may be considered the "outer reach" of the atmosphere. The results of 
these calculations are presented as the solid curves in figures 2 and 3* 

It is seen in figure 3 that for the high entrance speeds considered, 
the decelerations reach large values con?>ared to the acceleration of 
gravity. This suggests that the gravity term in equation (3) may be 
neglected without serio\isly eiffectlng the results When this term is 
neglected the equation of motion becomes 


_ ^ V ^ e"^ V2 

dt dy 2m 


( 8 ) 


®It is usual to neglect the gravity acceleration a priori (see e.g., 
refs. 1 and 14-.) 
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Integration, gives 


to V 


-py 

2pm 


+ const. 


or 


„ e“^ 

V = const. X e 

At the altitude of ifO miles it can readily he shown that the term 

-py 

® 

e spa 

is very nearly unity so that the velocity may he written 

_ ^PqA ^-py 

V = Vj. e 2pm (9) 

and 


d7 

dt 

e 


CdPoAVe^ -py - 

e e 

2mg 


^PPQ-^ -py 

pm 


( 10 ) 


where Vg is the entrance speed. 

By use of equations (9) and (lO) the vertical -decent speeds and 
decelerations for the 1-foot -diameter sphere previously considered have 
heen calcialated for the same entrance speeds. The results are shown as 
the dashed curves in figures 2 and 3. It is seen that these approximate 
calc^llations agree very well with those based on the more complete 
equation of motion (eq. (3))- 

ihe above finding is important, for it indicates that in the gen- 
eral case, wherein the body enters the atmosphere at high speed at 
angle 0g to the horizontal, the gravity term, provided is not 
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too small/ may "be neglected in equation (1) to yield 

^ sin 

dt^ 2m 

„ Cpp^A cos 
dt® “ 2m 


( 11 ) 


so that the fll^t path is essentially a straight line (l.e., 6 = 6g), 
and the resultant deceleration eqTjation hecomes 


dV 

dt 2m 


( 12 ) 


Now, again, if the density relation given by equation (2) is used and 
it is noted that 


dy 



sin 0 J 1 


dt 


V sin 6 


dV 
® dy 


equation (12) becomes 


dV _-3y, 

-==- = ————— e dy 
V 2m sin 


which can be integrated to yield 


V.Vj;e 

and the deceleration is then 

dV o e"^ 

dt CdPcAVe -Py ^ 

— - — at — ■ ■ e e 

g 2mg 


(13) 


( 11 ^) 
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T3ae altlinicLe yx "Whlcti the maximum deceleration occurs is found, 
from this relation to he 


Yi 


1 

P 


CpPpA 
^ sin 


(15) 


If Yx is positive, the velocity Yx (from eg.s. ( 13 ) and (15)) at 
■which the Tnaxlmim deceleration occurs hecomes 

Vi - Vj; e“^ = 0.61 Ve (16) 

and the val'ue of the maximum deceleration is 




(17) 


If equations (l3) and (l^) are rewritten to mahe the altitude 
reference point yx rather than zero, then 


_ g-P(yi+Ay) 

2 Bm. sin ©w 
V = Ve e ^ ^ 


and 


dv CpppA -3(yx-tAy) 

dt ^d'^o-^-Te^ -p(yi+Ay) Pm sin ©E * 

— = e e ^ 

S Pmg 

respectively, -where ^ is the change in altitude from yi. Substitu- 
tion of eq.uation (15) into these expressions can readily he shown to 
give 


^Ay 

'^E 


-i- 




= F’ ( PAy) 


(18) 
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and 



e-^ *F”(BSSr) 


(19) 


Equations (l8) and (19) are generalized expressions for velocity and 
deceleration for "bodies of constant drag coefficient and^ together with 
equations (15) and (17) ^ can be used to determine the variation of 
these quantities with altitude for specific cases. The dependence of 
F’(0;^) and F”(P^iy) on is shown in figure !(•. 

The maximum deceleration and the velocity for maximum decelera- 
tion as given by eq'uations (17) and (l6) apply only if the altittide yi, 
given by equation (15)# is positive. Otherwise the maximum decelera- 
tion in flight occurs at sea level with the velocity (see eq. ( 13 )) 




V = Vq = Vg e (20) 

and has the value 




2 ■■ CpPpA ... 

CpPpA-VE g pm sin 6^ 
2mg 


( 21 ) 


Heating of the Body 


It was noted previously that-for practicable rocket missiles# it 
is vital that the wei^t of the missile be kept to a minimum. The 
total heat transferred to a missile from the air must be absorbed by 
some "coolant” material. Since this material has a maximum allowable 
temperature# it follows that it can accept only a given amount of heat 
per •unit weight. Hence# the total heat input to the missile must be 
kept at a minimum for minimum missile weight. 
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Often the coolant material is sln^Jly the shell of the missile and 
as such must provide the structural strength and rlgldit 7 for the 
missile as well. Ihe strength of the structure is dictated, in pasrt, 
hy the stresses induced hy temperature gradients within the shell. 

Since these temperature gradients are proportional to the time rate of 
heat input, the maximum time rate of heat input is important in missile 
design. The heating, of course, varies along the surface hut, since 
the shell transmits heat along as well as through itself, the strength 
of the structure as a whole may he determined hy the maximum value of 
the average heat -transfer rate over the siirface. This is simply the 
maxInruTn value of the time rate of heat input per unit area. On the 
other hand, the structural strength at local points on the surface may 
he determined primarily hy the local rate of heat input. Hence, the 
maximum time rate of heat input per unit area at the surface element 
where the heat transfer is greatest may also he of importance in design. 

If liquid cooling is employed, the maximum surface heat-transfer 
rates retain their significance hut, now, in the sense that they dic- 
tate such requirements as maximum coolant pumping rate, or perhaps 
shell porosity as well in the case of sweat cooling. Whichever the 
case, in the analysis to follow, these elements of the heating problem 
will he treated: 

1. The total heat input 

2. The maximum time rate of average heat input per tinit area 

3. The 7Tw.-!fi rninri time rate of local heat input per unit area 

Since it is the primary function of this report to study means 
available to the missile designer to minimize the heating problem, the 
analysis is simplified to facilitate comparison of the relative heating 
of one missile with respect to another - accurate determination of the 
absolute heating of individual missiles is not attempted. With this 
point in mind, the following assumptions, discussed in Appendix B, are 
made: 

1. Convective heat transfer predominates (i.e., radiation effects 

are negligible) . 

2. Effects of gaseous imperfections may he neglected. 

3. Shock -wave houndary-layer interaction may he neglected. 

Reynolds’ analogy is applicable. 

5. Prandtl number is unity. 

Total heat input. - The time rate of convective heat transfer from 
the air to any element of surface of the hody may he expressed hy the 
well-known relation 


S - Tw)j 


( 22 ) 



12 


NACA BM A53D28 


where 

S heat transferred per unit area, ft-lb/ft^ 

ft-lh 

h convective heat-transfer coefficient, — 

f t^ sec °R 

Tj. recovery temperature , 

Tv temperature of the wall, °R 
t time, sec 

and the subscript 2 denotes local conditions at any element of the 
STJrf ace dS . 

It~is convenient in pairt of this analysis to determine the heating 
as a function of altittide. To this end, noting that 


dt 


-<3y 
V sin 


we see that equation (22) may he written 


M = _ ~ ^23) 

^ V sin 

With the asstm^jtion that the Prandtl number Is unity, the recovery tem- 
perature is 

where 

M Mach number at the altitude y, dimensionless 

7 the ratio of specific heat at constant pressure to that at constant 
volume, Cp/Cv, dimensionless 

T static temperature at the altitude y, °R 
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so that 




Tv)z = T - 


7-1 




T 


It is seen that for large values of the Mach ■ number , which is the case 
of principal interest, the third term is large compared to reasonably 
al.lowahle values of T - Tv It will therefore he asstnned that T - 
is negligible® so that 

(Tr - Tw)j = ^ M^T (2lf) 


Moreover, since 


M^ T = 


(r-i)Cp 


equation (24) may be written 


Now the local heat-transfer coefficient h^ is, by Reynolds' analogy, 
for the assimied Prandtl number of unity 

l>2=|CfjCpjPlVi ( 26 ) 


where C. 


f is the local skin-friction coefficient based on conditions 


Ihus, since (T^. - Tv) 
is essentially constant over the entire surface S, the rate of total 


V,, etc.. Just outside the boundary layer. 


®It shoTild be noted that without this assumption, the heat-input deter- 
mination would be greatly complicated since the changing wall tem- 
perature wiih altitude would have to be considered to obtain the heat 
input (see e.g., ref. l) . For high-speed missiles which maintain 
high speed during descent, the assuuqjtion is obviously permissible. 
Even for high-speed missiles which finally decelerate to low speeds, 
the asstimption is generally still adequate since the total heat input 
is largely determined by the heat transfer during the high-speed 
portion of fli^t. 
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heat transfer vlth altitiide 1)60011163 from eqmtions ( 23 ) through ( 26 ) 


<iQ r dH 
dy "7 <3y 


sin % 


/ J =Pi ‘’l ■5S 


where Q is the heat transferred to the whole surface S- This equa- 
tion may be written 

dQ Cf ‘PV^ S 

^ (27 

k sin 0g 

wherein Cp is set equal to Cp and 
“h 


Cf ’ = i r Cf ds 

sj ^2 V p y V V/ 


The parameter Cf* is termed ”the equivalent friction coefficient,” 
and will be assumed constantT^ independent of altitude, again on the 
premise that relative rather than absolute heating is of interest. 
With equations (2) and (l3 ) f then, equation (27) is written 


CpPpA • 

*^f '-^o^E^ "fiy ^ sin-0p. ^ 
4 sin 0-ri 


Comparison of eqmtion ( 29 ) with equation (l4) shows that the 
altitvide rate of heat transfer is directly proportional to the 


"^This assumption would appear poor at first- glance since the Mach number 
and Reynolds number variations are so large. Analysis has indicated, 
however, that the effects of Mach n\miber and Reynolds number variation 
are nearly compensating. The variation in Cf' for typical conical 
missiles was found to be, at most, about 50 percent from the maximtim 
Cf ’ in the altitude range in which 80 percent of the heat is 
transferred. 
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deceleration, so that 


dQ/ay 



mg / ^f \ 

2 sin 6 g V Cj)A / 


(30) 


and therefore the TTw-ytTmim altitude rate of heat transfer occurs at the 
altitude yi (see eq.. ( 15 )) and is given by 



(31) 


It follows, of course, that the altittide rate of heat transfer varies 
with incremental change in altitiode from yi in the same mann er as 
deceleration, and thus (see eq.. ( 19 )) 


(dQ/dy)^y 

(dQ/dy)^ 


F”(e^3) 


(32) 


The total heat input to the body at impact follows from equation ( 29 ) 
(integrating over the limits 0 < y < «) and is 


Q = 





CpPpA 

pm sin 0JJ 


The impact velocity, (the velocity of body at y = O), is 


(33) 


CpPoA 


2 an sin e. 

Vo = Ve e 


so that equation ( 33 ) may be written in the alternative form 



m 
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MaxlTmuB time rate of average heat Input per milt area. - To deter- 
mine the time rate of average heat transfer per imlt area> equations ( 25 )^ 
( 26 ) f and ( 28 ) vith equation (22) may he shown to give 


dHcv 1 a 

i Cf’ PV® 
dt 4 ^ 


(35) 


which, together with equations (2) and ( 13 )# hecomes at altitude y 


dt 




Cf'f’oV|] _Qy- 2^al Bin 0g 
: e ™ e 


(36) 


5he maxl.mum time rate of average heat transfer per \mit area is found 
from this expression to he 


(« 


■(%=) 


6e \CjjA/ 


mVi? sin 9 


E 


(37) 


and It occurs at the altitude 


3 20m sin 


(38) 


where the velocity is 


Ve e*" 3 sr 0.72 Ve 


(39) 


As with altitude rate of heat transfer, it can he shown that 


(dHav/dt)^ 

(dHav/dt)^ 


m 


F"0^y) 



E 
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Equations (3T)> (38), and (39) apply If the altitude for Tnaxirmnn time 
rate of average heat transfer per unit area occurs above sea level. 

It y 2 .f equation (38), is negative, then this rate occurs at sea 
level and is, from equation ( 38 ), 


spm Bln 

(4.1) 


Ifaximum time rate of local heat input per unit area. - Ihe elemen- 
tal surface -which is subject to the greatest heat transfer per vmit 
area is, except in unxisual cases, the tip of the missile nose which 
first meets the air. It seems unlikely that a poin-ted nose will be of 
practical interest for hi^-speed missiles since not only is the loc^ 
heat-transfer rate exceedingly large in this case, but the capacity 
for heat retention is small. Th-us a truly pointed nose woTild bum 
away. Body shapes of in-berest for high-speed missiles would more prob- 
ably, then, be those with nose shapes having nearly hemispherical tips. 
The following analysis applies at sTJch tips. 

It is well known that for any truly blimt body, the bow shock -wave 
is detached and there exists a stagnation point at the nose. Consider 
conditions at this point and assume that the local radl-us of curvature 
of the body is a (see sketch). 

The bow shock -wave is normal to 
the stagnation streamline and 
converts the s-upersonic flow 
ahead of the shock to a law- 
subsonic speed flow at high 
static temperat-ure downstream 
of the shock. Thus, it is s-ug- 
gested that conditions near the 
s-tagnatlon point may be investl- 
ga-ted by treating the nose sec- 
tion as if it were a segment of 
a sphere in a s-obsonlc flow field. 

The heat -transfer rate per 
tmit area at the stagnation point 
is given by the relation 




Nur ^r(^w " '^r) 



a 
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where is the thermal conductivity of the gas at the recovery tem- 

peratiare (i.e.^ total temperatiire) and Nur is the N\u 3 selt number 
of the flow. If the flow is assumed to be laminar and Incompressible^® 
BfUy is given^ according to reference 5 # Tsy the relationship 

_! _g 

Nur » 0.934 Pr ® 

We retain the ass\m 5 >tion that the Prandtl number is unity, note that 
RBjj = PVcr/lij., and substitute equation (25) into equation (kS) to obtain 


(43) 

Now it is we3J- known that at the high temperatures of interest here, 
the coefficient of viscosity Pj. varies nearly as the square root of 
the absolute temperature and is given by the relation 

M-r = 2.31 X 


If this expression is combined with equation (25) (neglecting T-^)^ 
equation (43) may then be written® 


dHs 

dt 


6.8 X I0‘®y| V® 


m 


®The assumption of constant density certainly may invalidate this anal- 
ysis for any quantitative st\idy of the relatively "cold -wall” flows 
of interest here. For the purpose of stu(3ying relative heat transfer 
g it should, however, prove adequate. 

Had wal.1 conditions rather than recovery conditions been used in the 
development of equation (44), the relation 

= 1-1 X 

would have been obtained assuming a linear variation of viscosity with 
temperature (to be consistent with the assumption of a cool wall) . 


This relation would give somewhat higher heat-transfer rates per unit 
area than equation (44) at velocities greater than about 36 OO feet-^oer 
second. 
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wbicTa, when combined with eqmtions (2) and (13) ^ becomes 


^ = 6.8 X 10“®y5 V/ e“ 


■— s^m sin 0-g 


( 45 ) 


The maximum valvie of dHg/dt can readily be shown to be 






BeCpOA 


(46) 


which occurs at the altitiide 


3 P 


= 1 ( ^f°_) 

p \ Pm sin 9-^y 


(47) 


corresponding to the velocity 


Vg = Vg e e a? 0.85 V- 


E 


(46) 


The manner in which the heat-transfer rate per unit area at the stagna' 
tion point varies with incremental change in altitude from, yg can be 
shown to be 


(fflsM)^ 

(dHs/dt)^ 




(^) 


(49) 


The dependence of F'''(p^) on is shown in figure 4. 

Equation (46) applies only if is above sea level. If yg, 
from equation (4?) ^ is negative, then the maximum heat-transfer rate 
per \init area at the stagnation point occurs at sea level and is 


'i'fc 'mas ^ 4 


aCpPoA 



(50) 
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DISCUSSION 

Motion 


The motion study shovs some important features about the high-speed 
decent of missiles throvigh the atmosphere. The major assumptions of 
this analysis -were that the drag coefficient was constant and the den- 
sity varied exponentially with altitude. It^was found that- the decel- 
eration due to drag was generally large compared tO“i;he acceleration 
of gravity and, consequently, that the acceleration of gravity cotild be 
neglected in the differential eq\iations of motion. The flight path was 
then seen to be a straight line, the missile maintaining the flight- 
path angle it had at entry to the atmosphere. 


For missiles presently of more usual interest, the maximum decel- 
eration occurs at altitude. One of the most interesting features of 
the flight of such a missile is that the maximum deceleration is inde- 
pendent“of physical characteristics (such as mass, size, and drag coef- 
ficient- of the missile), being dependent only on the entry speed and 
flight-path angle (see eq. (17)). The missile speed at maximum decel- 
eration (eq. (16)) bears a fixed relation to the entrance speed (6l per- 
cent of entremce speed), while the corresponding altitude (eq. (15)) 
depends on the physical characteristics and the flight-path angle but 
not on the entrance speed. It is also notable that for a given incre- 
mental change in altitude from the altitude for maximum deceleration, 
the deceleration and speed bear fixed ratios to the maximum decelera- 
tion and the entry speed, respectively (see fig. i)- and eqs. (I9) 
and (18)), hence, the deceleration and speed variation with altitude 
can readily be determined. 

If the missile is very heavy, the calculated altitude for maximum 
deceleration (eq. (15)) maybe fictitious (i.e., this altitude is neg- 
ative) so the maximum deceleration in flight, which occTors j^lst before 
impact at sea level, is less than that calculated by equation (l?) and 
is dependent on the body characteristics as well as the entry speed 
and flight-path angle (see eq. ( 21 )). However, the variation of speed 
and deceleration with altitude from the fictitious altitude given by 
equation (15) can still be obtained from figure if. 


Heating 


Total heat input . - In the heating analysis, a number of simplify- 
ing assumptions were made which shovild limit its applicability to the 
determination of relative values of heating at hypersonic speeds. It 
is in this relative sense that the following discussion pertains. 
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In considering the total heat transferred, hy convection to a 
missile, it is evident from eqiaation (33) that the ccnxrse the designer 
sho-uld take to obtain the least heating is affected hy the value of 
the factor 


CpPpA 
pm sin 


(51) 


To illustrate, first consider the case of a "relatively heavy" 
missile for which this factor is smal 1 compared to unity (the term 
"relatively heavy" is used to denote that the denominator involving 
the mass is very large as compared to the numerator involving the drag 
per unit dynamic pressure, CjjA) . Then 

CpPpA 

1 _ e sin ©g 


is sTTial 1 compared to 1. If this function is e 2 q>anded in series and 
only the leading term retained, equation (33) becomes 


Cf’SPoVE^ 
4p sin ©g 


(52) 


For the relatively heavy missile, then, the least heat will be trans- 
ferred when Cf'S is a minimum - that is to say, when the total shear 
force acting on the body is a minimum. This result is as would be 
expected, if one notes that requirixig B<<1 is tantamount to requir- 
ing the missile to be so heavy that it is retarded only slightly by 
aerodynamic drag in its motion through the atmosphere. Hence, the heat 
input to the missile is simply proportional to the shear force. 

Now let us consider the case when B > > 1, or, in other words, 
when this missile is "relatively light.” In this event, 

CpPpA 

1 _ e I3ni sin ©g ^ ^ 


and equation (3l) can be approximated 



( 53 ) 
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For the relatively light missile, then, the least convective heating is 
obtained when Cf'S/CjjA is a minimum. !Ihis is at first glance a 
rather surprising result, for it indicates that the heating is reduced 
by increasing the total drag, provided the equivalent frictional drag 
is not increased proportionately £is fast. Physically, this anomoly is 
resolved if the problem is viewed in the following way; The missile 
entering the atmosphere has the kinetic energy l/2 but, if 


CpPpA 
pm sin 



is small, then nearly all its entr^ce kinetic energy is lost, due to 
the action of aerodynamic forces, and must appear as heatirtg of both 
the atmosphere and the missile. The fraction of the total heat which 
is given to the missile is,^° from equation (33)# 


1 (Sill') 

2 \ CpA / 

Thus, by keeping this ratio a minimum, as much as possible of the 
energy is given to the atmosphere and the missile heating is therefore 
least. 


In order to ill-ustrate these considerations in greater detail, 
calculations have been made using the previously developed equations 
to determine the heat transferred by convection to a series of conical 
missiles. Two classes of missiles have been considered. Missiles in 
the first class were required to have a base area of 10 square feet. 
Missiles in the second class were required to have a volume of l6 cubic 
feet. Gross weights of 0, 1,000, 5,000, 10,000, and w pounds have 
been assumed, and the entrance angle, has been taken as 30 ° of arc 
in all cases. Missile heating, up to the time of Impact, has been 
calculated as a function of cone angle for entrance speeds of 10,000, 
20,000, and 30,000 feet per second. In these calculations the pressure 
drag coefficient was taken as constant for a particular cone at the 
value corresponding to the entrance Mach number (a valvie of T^, = 500° R 
was assumed throu^out) . These coefficients were determined from ref- 
erence 6 for cone angles of 10° and greater. For cone angles less 
than 10°, reference 7 was employed to determine these coefficients 


^°Note that even if all the drag is frictional drag, only half the heat 
is transferred to the body. The other half is contained in the 
bo^mdary la 3 rer and is left in the air in the body wake. 
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("base drag was neglected in al i cases) . The total drag coefficient 
was taken as the stim of the pressiire drag coefficient pltis the skin- 
friction coefficient, the latter coefficient being taken at its value 
for maximum, total heat-input rate with altitude. The houndary layer 
was assumed to he wholly turbulent since the Reynolds number, based on 
length of run along the sturface of a cone and local conditions just 
outside the boundary layer, was always greater than about 6 X 10® and, 
in fact, was of the order of billions for the more slender cones. 
Turbulent -boundary-layer data were obtained from references 8 and 9^ 
and Sutherland' s law for the variation of viscosity with temperature 
was used in obtaining ”eq.ul valent flat-plate” heat-transfer coefficients. 

Missile heating calculated in this manner for the fixed-base -area 
and fixed-volume cones is presented in figures 5 and 6, respectively. 
Curves for missiles having densities greater than steel are considered 
Improbable and are shown as dashed lines. It is clear that for both 
classes of bodies, when the missile is relatively heavy, the optimum 
solution is obtained by making C^'S as small as possible (small cone 
angle case) and this optimum is accentuated with increase in speed. On 
the other hand, when the missile is relatively light, reduced heating 
is obtained by making Cf'S/CjjA as small as possible (the large cone 
angle case) . It is noted also that, in general, the advantage of 
reduced heating of the relatively light, blunt cones is more pronounced 
in the fixed-base -area case than in the fixed-volume case. 

Maximum time rate of average heat input per unit area .- It was 
previously noted that. the maximum time rate of average heat input per 
tmit area may be of serious importance in determining the structural 
integrity of missiles entering the atmosphere at high speeds. n In 
order to illustrate this fact, consider the case of a missile having a 
shell made of solid material and assume that the rate of heat transfer 
per unit area does not vary rapidly from one surface element to the 
next. Then the rate of transfer of heat along the shell will be small 
courpared with the rate of transfer throu^ the shell. The shell stress 
due to heat transfer is that resulting from the tendency toward differ- 
ential expansion throu^ the shell and it is proportional to dT^/di] 
where T^^ is the temperature at any point t] within the shell and 
is measured perpendiciilar from the shell surface. We define k^ as 
the thermal conductivity of the shell material; then the rate at which 
heat transfers throu^ the shell per unit area is k^(dT.f|/dT]) and this 
must, at -r] = 0, equal the rate of heat input per unit surface area. 

For the missile considered as a whole, the maximum value of the average 
thermal stress in the shell is a measure of the over -all structural 


^^This is the common case when the shell material acts as structural 
support and must also transport or absorb the heat. 
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integrity and the maximum val^le of this stress will occur at the sur- 
face when 


_ 1 ^ ^ 
dt S \ dt y 

is a maximum. 

The course the designer should take to minimize the thermal stress 
for the missile as a whole is dependent, as for the case of total heat 
input, upon -vrtiether the missile is relatively heavy or light. For the 
relatively heavy missile the value of B, given hy equation (5l)^ is 
small compared to unity. The maximum value of the average therm^ 
stress in this case is proportional to (see eq. (4l)) 


. dt Jq ^ 




and, hence, the least average thermal stress is obtained by making Cf ' 
a minimum. On the other hand, for the relatively lig^t missile the 
maximum value of the average thermal stress is proportional to (see 
eq. (37)) 



(55) 


and, hence, the least average thermal stress occurs when Cf’/^DA 
is a minimum. 


In order to illiostrate these considerations in greater detail, 
the maximum values of the time rate of average heat input per unit area 
have been calculated for the const^t-base-area and the constant -volume 
cones previously discussed in the section on total heat- input. These 
values were determined in much the same manner as those of total heat 
input, with the exception that Cf* was evaluated at y^ (rather 
than yi), given by equation ( 38 ) when it applies, and otherwise at 
y = 0. The results are shown in figures 7 and 8. It is seen that 
the maximum values of average thermal stress are reduced for both the 
slender cones and blvint cones as con^jared to the relatively large 
values of this stress experienced by cones of Intermediate slenderness. 


Maximum time rate of local heat input per unit" area. - Perhaps even 
more important than the maximum value of the average shell stress is 
the maximum stress that occurs in the shell at the surface element of 


.E 
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■fche missile nose,^ ■wiiere tlie local liea'b-'transfer rate is proTaatly the 
greatest, for, in general, this latter stress is many times larger. 

In fact, this rate of local heat input can he so large as to promote 
tei^erature gradients throxi^ the shell that are intolerable even with 
the most highly conductive materials (copper, silver, etc.).^® Hhus 
some additional means of cooling, such as sweat cooling, may, in any 
case, be req.uired in this region. 

It was stated previously that pointed-nose bodies are undesirable 
due, in part, to the fact that the local heat -transfer rate per imit 
area at the tip is excessive. The validity of this statement is 
demonstrated by the results of the analysis. It is clear (see eq.. (44)) 
that since the local transfer rate varies inversely with the sqtiare 
root of the tip radius, not only should pointed bodies be avoided, but 
the rounded nose should have as large a radivis as possible. Rie ques- 
tion then arisesj if the nose radius is arbitrarily fixed, what course 
is available to the missile designer to minimize the problem of local 
heating at the stagnation point? From both equations (46) and (50), 
it is seen that for an arbitrary nose radius, if the mass, entry speed, 
and flight -path angle are fixed, then the only way to reduce the stag- 
nation rate of heat input per xmit area is to increase the product CpA. 
In fact, a relative stagnation-point heat-transfer rate per unit area, 
may be expressed In terms of B (see eq. (51))^ if it is defined as the 
ratio of the maximum stagnation-point heat -transfer rate per unit area 
for a given missile to the maximum rate the same missile would experi- 
ence if it were infinitely heavy. For the infinitely heavy missile, 
the maximum rate occurs at sea level and is (see eq. (50)) 



so that from equation (50) 


" 2$m slh e„ - 
+ = e ® = e ^ 


(56) 


In this report we are concerned only with bodies. If wings or sta- 
bilizers are used, their leading edges are similarly svirface elements 
which experience intense heat transfer. The heating problem with 
wings and stabilizers is, in fact, so serious at very high speeds 


that their use as lifting surfaces appears, at present, inadvisable. 
^®See reference 1 for further discussion. 
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if the given mlssiXe £Q.so attains its wAYiTmim rate at sea level 
(i.e., y-g = 0; eq,. (1}-T))j T(herea.s 


-/ 


pm sin 


SeCj^PgA ^ 2eB 


( 57 ) 


if the given missile attains its naxinrum rate ahove sea level (eq.. (46), 
y-g positive) . The variation of i|f with 1 /b is shown in f igi^e 9. 
Clearly, the high presstcre drag shape has the advantage over the slender 
shape in this respect. 

In order to illustrate these considerations in greater detail, 
again consider the constant -hase -area and constant -volume cones dis- 
cussed earlier. Assume the pointed tips of all the cones are replaced 
hy spherical tips of the same radius a. The relative effect of vary- 
ing the cone angle on the stagnation-point heating can then he assessed 
hy determining the variation of the product 



This product has heen calculated for the various cones, assuming Cjj 
to he una.ffected hy the addition of the hemispherical tip (the tip 
radius may he arbitrarily small), and the results are shown in figures 10 
and 11. It is seen again that the missiles having large cone angle 
(high drag coefficient) sire considerably superior. 


DESIGN C0NSIDERA.TI0NS AND CONCLUDING REMARKS 


In the foregoing analysis sind discussion, two aspects of the heat- 
ing problem for missiles entering the atmosphere were treated. The 
first concerned the total heat absorbed hy the missile and wsis related 
to the coolant required to prevent its disintegration. It was found 
that“if a missile were relatively light, the least required weight of 
coolant (and hence of missile) is obtained with a shape having a high 
pressure drag coefficient, that is to say, a blunt shape. On the other 
hand, it was found that if the missile were relatively heavy the least 
required weight of coolant, and hence of missile, is obtained with a 
shape having a low skin-friction drag coefficient, that is to say, a 
long slender shape. 

The second aspect of the heating problem treated was concerned 
with the rate of heat input, partlctaarly with regard to thermal shell 
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stresses resulting therefrom. It -was seen that the maximum average 
heat-input rate and, hence, maximum average thermal stress coxild he 
decreased hy using either a blunt or a slender missile, while missiles 
of Intermediate slenderness were definitely to be avoided in this con- 
nection. The region of highest local heat-transfer rate and, hence, 
probably greatest theimal stress was reasoned to be located at the 
forward tip of the missile in most cases. This was assumed to be 
the case and it was found that the magnitude of this stress was reduced 
by employing a shape having the largest permissible tip radius and 
over-all drag coefficient; that is to say, the blunt, hl^ drag shape 
always appears to have the advantage in this respect. 

Biese results provide us with rather crude, but useful, bases for 
determining shapes of missiles entering the atmosphere which have mini- 
mized heat-transfer problems. If the over-all design considerations 
of payload, booster, et al, dictate that the re-entry missile be rela- 
tively heavy in the sense of this report, then it may be most desirable 
to make this missile long and slender, especially if the entry speed is 
very high (say 20,000 ft/sec or greater). Perhaps the slender conical 
shape is appropriate for such a missile. It seems clear, too, that the 
tip of this missile should be given the largest practicable nose radius 
in order to minimize the maximum local heat-transfer rate and hence 
maximum local shell stress problem. Even then it may be necessary to 
employ additional means to minimize the heat -transfer rate and, hence, 
thermal stress- encountered in this region (e.g., by sweat cooling). 

Let us now consider the case where the over-all design conditions 
dictate that the re-entry missile be relatively light in the sense of 
this report. This case is believed to be of more immediate ln^jortance 
than the one just considered since the lower sizes, wei^ts, and 
entrance speeds -to which it applies are more nearly in line -with those 
presently of in-terest. The relatively light re-entry missile will 
therefore be treated at grea-fcer length. 

A shape which should -warrant at-tention for svich missile applica- 
tion is the sphere, for it has the following ad-vantages; 

1. It is a high drag shape and the frictional drag is only a few 

percent of the to-tal drag. 

2. It has the maximum volume for a given surface area. 

The continuo-usly curved surface is inherently stiff and strong. 

The large s-tagnation-point radi-us significantly assists in 
reducing the maximum thermal stress in the shell. 


3- 

k. 



28 


KACA EM A53D28 


5- Aerocfynamic “forces are not eeneltive to attitude and, hence, 
a sphere may need no stahilizing surfaces. 

6. Because of this Insensitivity to attitude, a sphere may pur- 
posely he rotated slowly, and perhaps even randomly^'* 
during flight. In order to svdjject all surface elements to 
about the same amount of heating and thereby approach uni- 
form shell heating. 

On the other hand, the sphere. In counnon with other very hl^ drag 
shapes may be unacceptable If: 

1. Ihe low tenninal speed permits effective countermeas-ures. 

2. The lower average speed of descent Increases the wind drift 

error at the target. 

3. The magnitude of the maximum deceleration Is greater than can 

be allowed. 

The first two of these disadvantages of the sphere ml^t be minimized 
by protruding a flow-separatlon-lnduclng spike from the front of the 
sphere to reduce the drag coefficient to roughly half ( see ref . ll) . 
Stabilization would now be required but only to the extent required to 
counterbalance the moment produced by the spike. Special provision 
wo^LLd have to be made for cooling the spike. 

These possible disadvantages of very hl^ drag shapes may also be 
alleviated by another means, namely, \islng variable geometry arrange- 
ments. For example, an arreingement which suggests Itself is a roTjnd- 
nosed shape with conical afterbody of low apex angle employing an 
extensible skirt at the base, as shown In figure 12. With the skirt 
flared, the advantages of high drag are obtained during the entry 
phase of fli^t. As the air density increases with decreasing altitude, 
the skirt flare is decreased to vary the drag so as to produce the 
desired deceleration and speed history. If the deceleration is speci- 
fied In the equation of motion (see motion analysis), the required 
variation of drag coefficient with altitude can be calculated and, in 
turn, the heating characteristics can be obtained. 


^■^ote that if rotation Is permitted, slow, random motion may be 
required in order to prevent Magnus forces from causing deviation 
of the flight path from the teirget. It should also be noted that 
at subsonic and low supersonic speeds gun-fired spheres, presumably 
not rotating, have shown rather large lateral motions in flight (see 
ref. lO) . It is not known whether such behavior occurs at high 
supersonic speeds. 
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The examples considered, of course, are Included only to demon- 
strate some of the means the designer has at hand to control and 
diminish the aerodynamic heating problem. For si ngilicity, this problem 
has been treated, for the most part, in a relative rather than absolute 
fashion. In any final design, there is, of course, no substitute for 
step-by-step or other more acctirate calculation of both the motion and 
aerodynamic heating of a missile. 

Even from a qualitative point of viev, a further word of caution 
must be given concerning the analysis of this paper. In particular, 
throu^out, we have neglected effects of gaseous imperfections (such 
as dissociation) and shock-wave boundary-layer interaction on convec- 
tive heat transfer to a missile, and of radiative heat transfer to or 
from the missile. One would not anticipate that these phenomena would 
significantly alter the conclusions reached on the relative merits of 
slender and blunt shapes from the standpoint of heat transfer at 
entrance speeds at least up to about 10,000 feet per second. It can- 
not tacitly be assinned, however, that this will be the case at hi^er 
entrance speeds (see Appendix B). Accurate conclxisions regarding the 
dependence of heat transfer on shape for missiles entering the atmos- 
phere at extremely hi^ supersonic speeds must await the availability 
of more reliable data on the static and dynamic properties of air at 
the hl^ temperatures and pressures that will be encountered. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Apr. 28, 1953 
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APPENDIX A 
SYMBOLS 


A reference area for drag evaluation, ft^ 

B body factor, dimensionless 

(See eq. (51)*) 

Cjj drag coefficient, dimensionless 

Cf skin-friction coefficient based on conditions just~outside 

the boundary layer, dimensionless 


Cf* 


Cp 


equivalent skin-friction coefficient, dimensionless 
(See eq. (28) .) 

ft -lb 

specific heat at constant pressure, 

slug 


Cv 


specific heat at constant volume. 


ft -lb 
slug °R 


F' ,F' ' ,F' ' ' 
g 


functions of dimensionless 

(See eqs. (18), (19) > and (^9)*) 

acceleration due to force of gravity ^taken as 32.2 



h 


convective heat-transfer coefficient. 


ft -lb 
ft^ sec °R 


H 

k 

m 

M 

Nu 


heat transferred per unit area. 


ft-lb 

ft^ 


thermal conductivity, ' ■ ' — 

sec ft^ (°R/ft) 

mass, slugs 

Mach number, dimensionless 
Nusselt number, dimensionless 
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Pr 

Q 

Re 

S 

T 

t 

V 

z 

P 

7 

A 

0 

P 

a 

t 


0 

1 


Prandtl n-uiriber, dimensionless 
total "heat transferred, ft-lb 
Resmolds nxcmber, dimensionless 
surface area, ft^ 

temperature (ambient temperatiare of air at altitude y 
unless otherwise specified) , °R 


time, sec 
velocity. 


ft 

sec 


horizontal and vertical distance from impact point, ft 

ft^ 

variable of Integration, 

sec^ 

constant in density - altitude relation, ft“^ 

(See eq. (2) . ) 

ratio of specific heat at constant pressure to specific 
heat at constant volume, Cp/c^, dimensionless 

increment 

distance within the shell measured normal to shell surface, ft 

angle of flight path with respect to horizontal, deg 

slugs 


coefficient of absolute viscosity. 


ft sec 


air density, 
radius, ft 


slug 

ft® 


relative heat-transfer factor, dimensionless 
(See eqs. (56) and (5?)-) 


Subscripts 

conditions at sea level (y = o) 
conditions at altitude yi(eq. (15)) 
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2 

3 

E 

1 

r 

s 

w 

n 


conditions at altit\ade y2(eq. (38)) 

conditions at altitude y„(eq,. (^T)) 

conditions at entrance to earth's atmosphere 

local conditions 

recovery conditions 

stagnation conditions 

wall conditions 

conditions within the shell of the missile 
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APPENDIX B 

SIMPLrFYTNG ASSUMPTIONS IN IHE CALCULATION OF 
AERODYNAMIC HEATING 


As noted In the main hody of the report, the heating analysis is 
simplified by making the following assumptions: 

1. Convective heat transfer is of foremost importance; that is, 

radiative effects may be neglected. 

2. Effects of gaseous imperfections, in particular dissociation, 

may be neglected. 

3. Effects of shock -wave boundary-layer interaction may be 

neglected. 

4. Reynolds' analogy is applicable. 

5- Prandtl number is unity. 

The restrictions imposed by these assumptions will now be considered 
in some detail. 

In assumption 1, two simplifications are involved; namely, (l) radia- 
tion from the surface of the body is neglected, and (2) radiation to the 
body from the high-te^perat^Jre disturbed air between the shock wave and 
the surface is neglected. The first simplification may be Justified on 
the premise that the maximum allowable surface temperature will be about 
the same for one body as compared with another, irrespective of shape, 
and, consequently, radiation away from the surface will be approximately 
the same. Hence, neglecting this form of heat transfer should not 
appreciably change the relative heating which is of principal interest 
in this paper. 

The second simplification of ignoring radiative heat transfer from 
the disturbed air to the body is not so easily treated. At ordinary 
flight speeds this form of heat transfer is negligible since it is well 
established that at temperatures not too different from ambient tem- 
perature, air is both a poor radiator and a poor absorber. At the 
flight speeds of interest, temperatiores in the tens of thousands of 
degrees Fahrenheit may be easily obtained in the disturbed air flow, 
especially about the heavier blunt bodies. At these temperatiures it 
does not follow, a priori, that air is a poor radiator. Data on the 
properties of air at these temperatures are indeed meager. Hence, it 
is clear that calculations of radiative heat transfer from air under 
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these conditions must, at "best, he qtialitative . Nevertheless, several 
such calculations have been made, assuming for lack of better informa- 
tion that air behaves as a grey body radlator_and that Vein’ s law may 
be used to relate the wave length at which the maximum amount of radia- 
tion is emitted to the temperatiare of the air (this assun 5 )tion, in 
effect, enables low-temperature data on the emissivity of air to be 
used in calculating radiation at hi^ temperatures) . In these calcu- 
lations effects of dissociation in reducing the temperature of the 
disturbed air have also been neglected and hence from this stsindpoint, 
at least, conservative (i.e., too high) estimates of- radiative heat 
transfer should evolve. The results of these ceilculations indicate 
the following: (l) Radiative heat transfer from the disturbed air to 

the body is of negligible importance compared to convective heat trans- 
fer at entrance speeds in the nei^borhood of, or less than, 10,000 feet 
per second; (2) Radiative heat transfer, in the case of relatively mas- 
sive hl\mt bodies, may have to be considered in heat-transfer calcula- 
tions at entrance speeds in the neighborhood of 20,000 feet per second; 
(3) Radiative heat transfer, in the case of relatively massive blunt 
bodies, may be of comparable importance to convective heat- transfer at 
entrance speeds in the neighborhood of 30,000 feet per second. From 
these resTilts, we conclude, then, that the neglect of radiative heat 
transfer from the disturbed air to the body is probably permissible for 
all except, perhaps, very blunt and heavy shapes at entrance speeds -up 
to 20,000 feet per second. However, this simplification may not be 
permissible, especially in the case of heavy blunt bodies entering the 
atmosphere at speeds in the neighborhood of, or greater than, 30^000 
feet per second. 

In assumption 2, the neglect of effects of gaseous imperfections, 
partlculeirly dissociation, on convective heat transfer would appear to 
he penaissible at entrance speeds up to and in the nei^borhood of 

10.000 feet per second, since at such speeds the temperatures of the 
disturbed air are not hi^ eno^lgh for these imperfections to become 
significantly manifest. On the other hand, as the entrance speeds 
approach 20,000 feet per second, temperatures of the disturbed air may 
easily exceed 10,000° Rankine, in which case appreciable dissociation 
may be anticipated, inside the boundary layer for all bodies, and 
inside and outside the boundary layer in the case of blunt bodies. The 
magnitude of these effects is at present in some doubt (see, e.g., the 
results of refs. 12 and 13- ) Hence, for the present, the neglect““of 
effects of gaseous Inperfections on convective heat transfer is not 
demonstrably permissible at entrance speeds in the neighborhood of 

20.000 feet per second or greater. 

In assumption 3^ it has been shown by Lees and Probsteln (ref. l4) , 
and more recently by Li and Nagamatsu (ref. 15)# that shock -wave 
boiindary-layer interaction may significantly increase laminar skin- 
friction coefficients on a flat plate at zero incidence and Mach 
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numbers in excess of about 10. Lees and Protstein fo-und somewliat the 
opposite effect on heat-transfer rate in the case of west interaction. 

It is not now known how this phenomenon depends upon hody shape or tj^pe 
of houndary layer. However, it is reasonable to anticipate that there 
will be some effect, and certainly if the skin-friction coefficient is 
increased in order of magnitude at Mach numbers approaching 20, as 
indicated by the res\ilts of Li and Nagamatsu for strong interaction, 
then the phenomenon cannot be presuBned negligible. Hence, we conclude 
that from this standpoint, also, the convective heat-transfer calcula- 
tions of this report may be in error at entrance speeds of the order 
of 20,000 feet per second or greater. 

Ihe assTimption that Reynolds' analogy may be used to relate skin- 
friction and heat-transfer coefficient does not, especially in the 
light of recent work by Rubes in (ref. 16), seem out of line with the 
pmposes of this paper, at least at entrance speeds to and in the 
nei^borhqod of 10,000 feet per second. However, it does not follow, 
a priori, that this assumption remains valid at substantially hi^er 
entrance speeds, especia3J.y in view of the imperfect gas and shock-wave 
boundary-layer-interaction effects already discussed. 

The assunrption of Prandtl number equal to unity would also appear 
permissible for the analysis of relative heating of missiles at the 
lower entrance speeds considered here- However, in view of the ques- 
tionable effect (see again refs. 12 and 13) of dissociation on Prandtl 
n'umber, it is not clear that this asstnnption is strictly valid at the 
intermediate and higher entrance speeds treated in this report. 

From these considerations it is concluded that the simplifying 
assumptions made in the main heat-transfer analysis of this paper will 
not significantly Influence the results at entrance speeds in the 
neighborhood of or less than 10,000 feet per second. However, at 
entrance speeds in the neighborhood of and greater than 20,000 feet 
per second, these res-ults must be viewed with skepticism. More accurate 
calculations of heat transfer at these speeds must, among other things, 
await more accurate determinations of both the static and dynamic prop- 
erties of air under these circumstances. 
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Figure 2 -Variations of velocity with attitude for a i-foot diameter, solid iron 
sphere entering the earth's atmosphere vertically at velocities of tO,OOQ 
20,000, and 30,000 ft /sec. 
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Figure 3 —Variations of deceleration with altitude for a I -foot diameter, solid 
iron sphere entering the earth's atmosphere vertically at velocities of 10,000, 
20,000, and 30,000 ft /sec. 
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Figure S." Convective heat transferred at impact to comcrd missiies of same base area entering the earth^ 

atmosphere at an ao^ of 30“ to the horizontal and velocities of tOfiOO, 20,000 and 30.000 ft /sec 
(base area » iO sq ft). 
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Figure 5.— ContmecL 
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Figure 5r Concluded. 
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Figure 6 Convective heat transferred at impact to conical missiles of same volume entering 
the earth's atmosphere at an angle of 30* to the horizontal and velocities of 10,000, 
ZOfiOO, and 30,000 ft /sec (volume • 16.34 cu ft). 
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Figure 6- Concluded. 
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Figure 7.- Maximum average rate of convective heat transfer to conical missiles of the same base area 
entering the earth's atmosphere at an angle of 30° to the horizontal and velocities of 10,000, 
20,000, and 30,000 ft /sec (base area • tO sq ft). 
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Figure Moxinum overoge rote of convective heot fronsfier to con/cot tnissifes of the some vo/urne 
entering the earWs atmosphere at an angle of 30^ to the horizontal and velocities of 10,000, 
20,000, and 30,000 ft/sec (volume - 16.34 cu ftX 
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Figure fO.- Maximum rate of convective heat transfer to the stt^nation point of spherically tipped cones 
of the same base area entering the earth's atmosphere at an angle of 30" to the horizontal and 
velocities of 10,000, 20,000, and 30,000 ft /sec ( base area = tO sq ft). 
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Figure 10.— Concluded. 
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Figure II. - Maximum rate of convective heat transfer to the stagnation point of spherically tipped cones 
of the same volume entering the earth’s atmosphere at an angle of 30” to the horizontal and 
velocities of 10,000, 20,000, and 30,000 ft /sec ( volume * 16.34 cu ft). 
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Figure //.- Continued. 
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